A Two-dimensional Map of the Color Excess in NGC 3603 
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ABSTRACT 

Using archival HST/WFG3 images centered on the young HD 97950 star cluster in the giant 
H II region NGC 3603, we computed the pixel-to-pixel distribution of the color excess, E(B — V) e , 
of the gas associated with this cluster from its Ha/Pa/3 flux ratio. At the assumed distance of 
6.9 kpc, the resulting median color excess within 1 pc from the cluster center is E(B — V) g = 
1.51 ± 0.04 mag. Outside the cluster (at r > 1 pc), the color excess is seen to increase with 
cluster-centric distance towards both North and South, reaching a value of about 2.2 mag at 
r = 2 pc from the cluster center. The radial dependence of E(B — V) s westward of the cluster 
appears rather flat at about 1.55 mag over the distance range 1.2 pc < r < 3 pc. In the eastern 
direction, E(B — V) g steadily increases from 1.5 mag at r = 1 pc to 1.7 mag at r = 2 pc, and stays 
nearly constant at 1.7 mag for 2 pc < r < 3 pc. The different radial profiles and the pixel-to-pixel 
variations of E(B — V) s clearly indicate the presence of significant differential reddening across 
the 4.9 pc x4.3 pc area centered on the HD 97950 star cluster. We interpret the variations of 
E(B — V) g as the result of stellar radiation and stellar winds interacting with an inhomogeneous 
dusty local interstellar medium (ISM) whose density varies spatially. From the E(B — V) s values 
measured along the rims of the prominent pillars MM1 and MM2 in the southwest and southeast 
of the HD 97950 cluster we estimate an H 2 column density of \og w (Nn 2 ) — 21.7 and extrapolate 
it to log 10 (A f H 2 ) = 23 in the pillars' interior. We find the pillars to be closer to us than the central 
ionizing cluster and suggest that star formation may be occurring in the pillar heads. 
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1. Introduction 

The young massive HD 97950 star cluster in 
the giant H II region NGC 3603 is located in the 
Sagittarius-Carina arm of the Milky Way. This 
cluster is one of the most massive young star 
clusters in the Galaxy (~ 10 4 M Q , Harayama et 
al. 2008). The cluster hosts ten times more OB 
stars than the Orion Nebula Cluster (ONC), in- 
cluding the two most massive binaries currently 
known in the Galaxy (Schnurr et al. 2008). The 
HD 97950 cluster displays pronounced mass seg- 
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regation (e.g., Sung & Bessell 2004; Harayama 
et al. 2008). Earlier studies trying to age-date 
the HD 97950 cluster found the massive stars on 
the upper main sequence (MS) to be very young 
(e.g., 1-2 Myr according to the spectroscopic 
study by Melena et al. 2008), while the pre-main- 
sequence stars (PMS) were found to show a larger 
age spread of 2 - 3 Myr (e.g., Eisenhauer et al. 
1998; Grebel 2004, 2005; Harayama et al. 2008; 
Pang et al. 2010). The PMS stars in the outer clus- 
ter regions and in the surroundings of the cluster 
may, in part, be older, but different authors arrive 
at different age ranges (e.g., 4-5 Myr; Sung & 
Bessell 2004; Rochau et al. 2010) or up to 10 Myr 
(Beccari et al. 2010). Also the evolved supergiants 
around the HD 97950 cluster support the idea of 
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a large age spread or alternatively sequential or 
multiple episodes of star formation (e.g., Moffat et 
al. 1983; Tapia et al. 2001; Crowther et al. 2008; 
Melena et al. 2008). On the other hand, the size 
of the wind-blown bubble around the cluster sug- 
gests an age much lower than 2.5 Myr when taking 
the kinetic energy input of the massive stars into 
account (see discussion in Drissen et al. 1995). At 
present, the age(s) and the star formation history 
of the HD 97950 star cluster and its surroundings 
in NGC 3603 are uncertain. One of the uncertain- 
ties in stellar age dating is introduced by the high 
and variable reddening in NGC 3603. 

The HD 97950 star cluster lies in a wind-blown 
cavity north of a large molecular cloud in the giant 
Hn region NGC 3603 (e.g., Clayton 1986; Melnick 
1989). The gaseous surroundings of the cluster 
show a complex and variable velocity and den- 
sity structure (e.g., Clayton 1990; Drissen et al. 
1995). Since the nebular density varies spatially, 
one would expect that also the dust extinction 
changes with position across the cluster area as 
denser clouds can shield dust from stellar radiation 
better than less dense clouds. The resulting vari- 
able reddening is known as differential reddening 
and was shown to be present in NGC 3603 (e.g., 
Sagar et al. 2001). When this effect is not taken 
properly into account, it may make stars appear 
redder and hence older in the color-magnitude di- 
agram (CMD) depending on their spatial position 
with respect to the line of sight. This may intro- 
duce a considerable uncertainty in the estimation 
of the stellar ages. 

Differential reddening could be one of the 
causes why the age spread among PMS stars in 
the HD 97950 star cluster appears to be as large 
as up to 10 Myr. In fact, as shown by Beccari et 
al. (2010, their Figure 7), PMS stars are spatially 
more widely distributed than the MS stars and 
reside in areas with different reddening. Sung & 
Bessel (2004) find that the color excess of the stars 
within the HD 97950 cluster core (at radii r < 0.7 
pc; Harayama et al. 2008) is E(B - V) s = 1.25 
mag and rapidly increases to 2.1 mag in the outer 
regions (r ~ 12 pc). 

We take advantage of the publicly available im- 
ages of the HD 97950 cluster taken with the Wide 
Field Camera 3 (WFC3) aboard the Hubble Space 
Telescope (HST) through narrow-band filters cen- 
tered on the Ha and Pa/3 lines and through broad- 



band filters sampling the continuum emission, in 
order to extend the work of Sung & Bessel (2004). 
We compute a pixel-to-pixel map of the color ex- 
cess, E(B-V) g , across a 2.43' x 2.14' (4.9x4.3pc) 
field centered on the HD 97950 cluster. Since only 
some of the gas contributing to the reddening is 
located in front of the majority of the stars, the 
reddening of the gas that we derive is an upper 
limit to the true reddening experienced by the 
stars. Thus the two-dimensional distribution of 
the color excess associated with the gas is the first 
step towards correcting individual MS and PMS 
stars for reddening in the WFC3 field of view. 
Multi-band stellar photometry will be needed in 
order to establish the conversion between the red- 
dening of the gas and that of the stars. This will 
ultimately allow us to constrain more tightly the 
age spread of PMS stars and the recent star for- 
mation history of the HD 97950 star cluster and 
of its immediate surroundings (Pang et al. 2011, 
in preparation). Our present paper focuses on the 
reddening map. The data and their reduction are 
described in Section 2, while the two-dimensional 
map of the color excess is presented and discussed 
in Section 3. Conclusions and a summary are pre- 
sented in Section 4. 

2. Observations and Data Reduction 

Multi- wavelength imaging of the HD 97950 star 
cluster was carried out in 2010 with fi5T/WFC3 
(proposal ID: 11360, PI: Robert O'Connell). For 
our analysis we use the optical images taken in 
the F555W (~ V), F656N (Ha), and F814W (~ 
I) filters with the WFC3 UVIS detector (0.04 
arcsec pixel -1 ) and near-infrared images in the 
F127M (continuum in J), F128N {Pa/3), and 
F139M (continuum in J) filters obtained with the 
WFC3 IR detector (0.13 arcsec pixel" 1 ). The 
total exposure time in the different filters is sum- 
marized in Table 1. 

All images were reduced with the WFC3 
pipeline and the IRAP0 MULTIDRIZ ZLE task. 
We note that no narrow-band filter sampling the 
continuum emission close to the Ha or Pa/3 line 
was used during the observations. Thus we use the 
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broad-band V and I filters to estimate continuum 
in Ha and the medium-band J-continuum filters 
for Pa/3. 

2.1. Ha and Pa/3 Emission 

Radiation emitted at shorter wavelengths is 
absorbed by dust more effectively than by pho- 
tons of longer wavelengths. Given that the Ha 
line (restframe wavelength 6563 A) is emitted at 
a shorter wavelength than Pa/3 (restframe wave- 
length 12802 A), the interstellar dust along the 
line of sight decreases the Ha flux detected by the 
observer more than the Pa/3 flux emitted by the 
same source. Therefore, the observed HajVs.fi flux 
ratio is smaller than the theoretical one computed 
for the same conditions of electron density (N e ) 
and temperature (T e ) of the source in the absence 
of dust. 

Equation (1) (taken from Calzetti et al. 1996) 
illustrates the relation between the color excess, 
E(B — V) g of the interstellar gas, and the observed 
and theoretical Ha/ Pa/3 flux ratios of this same 
gas (i? bs arid i?; n t, respectively) under the as- 
sumption of a specific extinction law (represented 
by k(X)=A(X)/E(B - V)): 



E(B - V) 9 



-log(R obs /R mt ) 
0A[n(X Ha ) - n{X Pali )] 



mag (1) 



In order to apply Equation (1), we need to de- 
rive the Ha and Pa/3 emission fluxes per pixel from 
the available WFC3 images, and for this purpose 
we make use of IRAF standard routines. We first 
determine an average point-spread function (PSF) 
in each filter for 15-20 stars in common to all im- 
ages. Since the PSF is largest in the F139M im- 
age, we degrade all the other images by convolving 
them with a Gaussian function whose dispersion, 
(T 2 , is the difference between the PSF a 2 of the 
F139M image and that of the image in the filter 
in question. We also re-scale the optical images 
to the same pixel scale of the near-infrared ones, 
and align all images to the one taken through the 
F139M filter. We calibrate all images in flux by 
multiplying them by their respective filter PHOT- 
FLAM value as given in the WFC3 manual (Drcs- 
sel et al. 2010). We derive the continuum emis- 
sion at the Ha (Pa/3) wavelength by interpolating 
the flux in the F555W and F814W (F127M and 



F139M) images pixel by pixel with a simple first- 
order polynomial. Such an interpolation allows 
us to take into account the slope of the continuum 
emission of the stars and to better remove the stars 
from the final, pure line-emission images. 

The continuum emission derived by inter- 
polating the flux between F555W and F8I4W 
(F127M and F139M) is subtracted from the 
F656N (F128N) image, and the output is mul- 
tiplied by the width of the narrow-band filter. In 
this way we are able to construct the images of 
the pure Ha and Pa/3 emission as well as their 
observed i? b s flux ratio. We derive i?; nt = 17.546 
from Osterbrock (1989) under the assumption 
that T e = 10 000 K and N e = 100 cm" 3 . We 
adopt a normal extinction law with a ratio of to- 
tal to selective extinction of Ry = 3.1, and derive 
n(X Ha ) = 2.355 and n(X Pal3 ) = 0.7644 from Fitz- 
patrick's (1999) extinction law. Finally, we apply 
Equation (1) to the image of the observed i? Q b s 
flux ratio to construct the pixel-to-pixel map of 
E(B — V) g shown in Figure 1. In panel (a) of Fig- 
ure 2 we show the histogram of the E(B — V) g per 
pixel, normalized by the total number of pixels in 
the i? bs image. The total line-of-sight color excess 
is always larger than the foreground reddening of 
E(B -V) = 1.1 mag (Pandey et al. 2000). More 
than 90% of the pixels have E(B — V) g between 
1.6 and 2.2 mag. 

2.2. Uncertainty of E(B - V) g 

In order to estimate the accuracy of our derived 
E(B — V) g values, we assume that the electron 
counts in the images follow a Poisson distribution. 
Therefore we start with a Poisson error on the 
electron counts per pixel in the input images, and 
propagate these errors through the entire compu- 
tation of the continua, their subtraction from the 
F656N and F128N images, and the computation 
of the E(B — V) g image according to Equation (1). 

We present the histogram of the uncertain- 
ties <Je(b-v) s P er pixel in panel (b) of Fig- 
ure 2, normalized by the total number of pix- 
els in the i? bs image. The histogram peaks at 
(?e(b-v) = 0-1 mag, and the percentage of pix- 
els with (Te(b-v) > 0-2 mag is less than 10%. In 
panel (c) of Figure 2 we plot the average (?e(b~V) & 
per pixel in bins of E(B — V) g per pixel. As can be 
seen, (Je(b~v) s gradually decreases with increas- 
ing E(B - V) g from 0.4 mag at E(B - V) g = 1.2 
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mag to 0.1 mag at E(B — V) g — 2.4 mag. 

In addition to photon noise, there are also some 
systematic errors because of the following effects: 

(1) Broad-band and medium-band filters are con- 
taminated by emission lines. The Pa/3 emission 
falls in the F127M filter. Assuming that the flux 
density of the continuum is the same in both the 
F127M and F128N filters, we can use the flux ra- 
tio F127M/F128N (where both filters were multi- 
plied by their respective PHOTFLAM and band 
width) to derive the contribution of the Pa/3 line 
to the flux in F127M. We find this contribution 
to be about 13%. No hydrogen lines are found 
in the wavelength range of the F139M filter, only 
a few faint He emission lines, which should not 
significantly contribute to the flux in this filter. 
The F555W filter includes four strong emission 
lines: the H/3 and the [Om]A = 4959 - 5007 A 
line, where the filter response is <~ 27% and the 
Ha line at a response of ~ 3%. The F814W fil- 
ter contains the [S ill] and [Ar ill] emission lines 
at wavelengths where its response is lower than 
10%. Unfortunately, we cannot estimate the con- 
tribution of these lines to the flux detected in each 
filter, because the only available spectrum taken 
of NGC3603 in the range of 3000 - 10400 A by 
Garcia-Rojas et al. (2006) was acquired at a po- 
sition outside the WFC3 field of view. If we as- 
sume that the line contamination is about 10% in 
F555W and F814 W as it is for F127M, the color 
excess would then increase by 0.05 mag on aver- 
age. 

(2) The adopted extinction law also contributes 
to the systematic uncertainties. If we replace 
Fitzpatrick's (1999) extinction law with that of 
Cardelli et al. (1989), we obtain an E(B — V) g sys- 
tematically smaller by 0.1 mag (Ry — 3.1). When 
we vary Ry by ±0.5 in Cardelli et al.'s (1989) ex- 
tinction law, E(B — V) g decreases by ~0.3mag 
(Ry = 3.6) or increases by ^0.2 mag {Ry = 2.6). 

(3) The assumed electron temperature and den- 
sity of the gas may be affected by systematic er- 
rors. Garcia-Rojas et al. (2006) and Lcbouteillcr 
et al. (2008) obtained T e = 10 000 K and 
N e = 1000 cm" 3 for the NGC3603 giant Hn 
region. Since this electron density is not avail- 
able in Osterbrock (1989), we calculated i?; n t for 
T e = 10 000 K and N e = 10 000 cm~ 3 and derived 
a new pixel-to-pixel map of E(B — V) g . A fac- 
tor of 100 difference in N e results in an average 



difference of 0.0025 mag in E(B — V) g , which is 
well within the errors in E(B — V) g due to photon 
noise. 

3. A Two-Dimensional Map of the Color 

Excess 

3.1. Global Properties of E(B - V) g 

The pixel-to-pixel map of E(B — V) g obtained 
for T e = 10 000 K and N e = 100 cm" 3 is presented 
in Figure 1. The resulting E(B — V) g is integrated 
along the different lines of sight towards NGC 3603 
and includes the Galactic foreground reddening, 
which amounts to E(B — V) = 1.1 mag accord- 
ing to Pandey et al. (2000). The HD 97950 cluster 
core (r < 0.52 pc) where dozens of OB stars reside 
is masked out in order to avoid saturated bright 
stars in the F555W and F814 W filters. Since the 
PSFs of the different filters may be slightly differ- 
ent even after the convolution procedure to obtain 
the same PSF described in Section 2.1, the sub- 
traction of the continua from the Ha and Pa/3 im- 
ages can produce negative values around the stars 
in the pure emission images, and thus negative 
E(B — V) g in the color excess map. This is par- 
ticularly true for saturated stars and their spikes. 
Therefore, we smooth the color excess map by re- 
placing the negative E(B — V) g at the position of 
a star with the median value of E(B — V) g in an 
annulus around it. 

In Figure 1, E(B — V) g is seen to decrease by 
0.2 mag when going from the eastern or the west- 
ern edge of the field of view to the cluster center. 
Overall, the gas color excess in the east is 0.1 mag 
larger than in the west. However, in the north- 
south direction E(B — V) g is systematically larger 
by 0.4 - 0.5 mag than in the east- west direction. In 
general, the southern region has the largest color 
excess as may be expected since here we are mov- 
ing towards the densest regions of the giant molec- 
ular cloud in the cluster vicinity (e.g., Niirnberger 
et al. 2002). Our findings support earlier sugges- 
tions that the line-of-sight dust density is uneven 
across the HD 97950 cluster surroundings. 

As an alternative way of quantifying the dif- 
ferential reddening around the HD 97950 cluster, 
we show in Figure 3 the radial dependence of 
E(B — V) g from the cluster core towards the east, 
west, north, and south. We define a series of con- 
centric annuli of the same width (0.2 pc) and with 
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increasing distance from the cluster center, and we 
measure the median value and the 16 and 84 per- 
centile values of the E(B — V) g distribution within 
each annulus. The solid lines in Figure 3 trace the 
median E(B — V) g as a function of cluster-centric 
distance, while the long-dashed lines indicate the 
84 and 16 percentile values. Towards the west, 
the color excess increases by 0.1 mag from 1 pc 
to 1.2 pc and stays rather flat at about 1.55 mag 
over the distance range 1.2 pc < r < 3 pc. In 
the eastern direction, E(B — V) g steadily increases 
from 1.5 mag at r = 1 pc to 1.7 mag at r = 2 pc 
and then remains nearly constant at 1.7 mag for 
2 pc < r < 3 pc. The most dramatic changes in 
E(B — V) g are seen along the northern and south- 
ern directions, where E(B — V) g gets larger by 0.3 
to 0.5 mag as the distance from the cluster center 
increases. Specifically, the color excess increases 
from about 1.7 mag at r = 1 pc to 2.2 mag at 
a distance of 2 pc. The pronounced increase to- 
wards the south does not come as a surprise since 
the overall density of the giant molecular clouds 
in NGC3603 increases in this direction (see, e.g., 
Figure 1 in Brandner et al. 1997a). 

In the two-dimensional map of the color ex- 
cess (Figure 1), some specific sources are associ- 
ated with a locally higher E(B — V) g . This is the 
case for two of the "proplyd" -like (protoplanetary- 
disks-like) objects first detected by Brandner et 
al. (2000). Proplyd 1 is located east of the 
HD 97950 star cluster, while Proplyd 3 is in the 
northwest. Both are tadpole-shaped and rim- 
brightened sources with extended tails pointing 
away from the cluster core. 

Mid-infrared observations (Niirnberger & Stankc 
2003) do not reveal any point-like source that 
might be associated with the proplyds. These au- 
thors suggest that these objects are small dense 
clumps of gas and dust, which are being photoe- 
vaporated by the intense ionizing radiation of the 
massive stars in the HD 97950 cluster instead of 
being disks around young stars as proposed by 
Brandner et al. (2000). In particular, Proplyd 3 
shows extended faint emission at 11.9 /<m, which 
may be caused by carbon-rich dust grains (see 
Goebel et al. 1995). Assuming that the dust prop- 
erties are the same for these three proplyds, the 
higher color excess of Proplyd 3 (0.3 mag redder 
than Proplyd 1, see Table 2) would be consistent 
with a higher dust content (if both proplyds are at 



the same distance), and hence with the fact that 
Proplyd 3 is the only proplyd detected at 11.9 ^im 
(Niirnberger & Stanke 2003). 

To the north of the HD 97950 cluster, we find 
Sher 25, a blue supergiant with a circumstellar 
ring and an hourglass nebula (Brandner et al. 
1997a, 1997b). We have measured the color excess 
in the bright rims of the northeastern and south- 
western lobes of the hourglass nebula (Table 2; see 
also Figure 1) and found that the hourglass nebula 
is more strongly cxtinctcd in the northeast (with 
E(B-V) g = 2.0 mag versus E(B-V) g = 1.6 mag 
in the southwest). 

3.2. The Shell Structure 

Optical images of NGC 3603 show an extended 
cavity around the HD 97950 cluster. Clayton 
(1986) estimated the diameter of this shell-like 
structure to be approximately 2 kpc. Clayton an- 
alyzed the Dopplcr shifts of the Ha emission in the 
central region of NGC 3603, and found a complex 
velocity structure in the ionized hydrogen. He in- 
terpreted the main velocity components as indica- 
tive of two expanding shells. One of these is a 
faint shell in the east-west direction with a mean 
heliocentric velocity of (f)hci = +120 km s _1 , an 
expansion velocity of ~ 30 km s _1 , and a diam- 
eter of ~ 1 pc. The other one is a brighter shell 
in the north-south direction and centered on the 
HD 97950 cluster. For this second shell-like feature 
Clayton (1986) found (w)hoi = +50 km s _1 and 
<~ 1.5 pc as a lower limit for its diameter. (Note 
that Clayton assumed that NGC 3603 lies at a dis- 
tance of 6 kpc from us.) From Clayton's Figure 
4, we infer an expansion velocity of ~ 55 km s _1 
for this second shell, which Clayton identified with 
the above mentioned cavity. 

In the E(B — V) g map (Figure 1) we also find a 
cavity in correspondence with this expanding gas 
shell and with a radius of about 1.0 pc. The me- 
dian color excess within the cavity is E(B — V) g = 
1.51 ± 0.04 mag. Moreover, we find a shell struc- 
ture corresponding to a local enhancement of color 
excess with a mean E(B — V) g — 1.59 mag about 
1.2 pc west of the cluster center (see Figure 4). 
This shell structure extends for 1.3 pc along the 
north-south direction (Figure 1). Judging from its 
position, it seems to lie along the high ionization 
front of the gas shell structure found by Clayton 
(1986) in the north-south direction. 
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Recently, Lebouteiller et al. (2007) observed 
NGC 3603 with Spitzer in order to search for PAH 
(polycyclic aromatic hydrocarbon) and [S iv] emis- 
sion. They found that within 1.5 pc from the 
cluster center, the PAH and [S iv] emission is sup- 
pressed, possibly because the strong radiation field 
of the massive stars in the cluster core (about 60 
OB stars, see Melena et al. 2008) destroys dust 
grains, dissociates complex molecules and pushes 
the gas out. In doing so, these stars would thus re- 
duce the color excess within the area and produce 
a cavity with lower E(B — V) g in the reddening 
map. The western edge of this cavity appears to 
be delineated by the shell structure seen in the 
reddening map. 

3.3. The Molecular Pillars 

In the southern region of NGC 3603 we can see 
two prominent pillars to the southwest and south- 
east of the cluster (Figures 1 and 4). Following 
Niirnberger et al. (2002), we refer to them as MM1 
and MM2, respectively. Figure 4 is a composite 
color image showing the Ha emission (in blue), 
the Pa/3 emission (in green), and the color excess 
(in red). The heads of both MM1 and MM2 stand 
out because of their strong Ha and Pa/3 emission. 
They are gradually being photoevaporated by the 
cluster's massive OB stars as is also indicated by 
the shocked and ionized material in their heads 
(Niirnberger et al. 2002, 2003). 

According to Bcrtoldi (1989) and Bertoldi & 
McKee (1990), in clouds experiencing photoevap- 
oration a velocity gradient will eventually emerge, 
leading to a configuration where the head of a pil- 
lar moves more slowly than its tail. Applying this 
to the Eagle nebula (M16), Pound (1998) points 
out that the measured radial velocity gradients de- 
pend on our viewing angle and provide clues about 
the three-dimensional structure of the nebular fea- 
tures with respect to the ionizing stars. The tail 
of a pillar located in front of the ionizing stars will 
then have lower radial velocities than the pillar's 
head, since the tail is being blown toward the ob- 
server (while the opposite radial velocity gradient 
would be observed across a pillar located behind 
the ionizing OB stars). Both Niirnberger et al. 
(2002) and Rollig et al. (2010) found that the ra- 
dial velocity in the head of MM1 is larger than 
that of the tail. Thus in NGC 3603 we have a con- 
figuration where the pillar is located in front of 



the stars of the HD 97950 cluster as seen from our 
position, with some tilt with respect to the light 
of sight. 

The rims of MM1 and MM2 have a high E(B - 
V) g of ~ 1.8 mag, implying a large amount of 
dust. Along the pillar's main body, though, the 
color excess is lower than at the rim (<~ 1.5 mag) 
owing to "limb brightening" effects. We explain 
this phenomenon in a cartoon in Figure 5. The 
pillar is represented by a cylinder whose density 
decreases from the head to the tail (Mackey & 
Lim 2010). The radiation of the cluster's OB stars 
ionizes the outer layer of the pillar (A, B and C in 
Figure 5), while the interior is largely shielded and 
mostly neutral. Because of the long path traveled 
by the light emitted in A and B of the ionized layer 
(Figure 5), the color excess measured in the pillar 
rims is relatively high, yet a lower limit to the true 
E(B — V) g associated with the main body of the 
pillar. The light emitted in C (i.e., the front side 
of the pillar in Figure 5) travels a shorter path 
within the ionized layer, and, for this reason, the 
E(B — V) g measured in the front side of the pillar 
facing the observer is lower than that computed 
at the pillar rims. 

Also in the simulations of Gritschneder et al. 
(2010) the ionizing radiation emitted by a mas- 
sive star can enhance the initial turbulent density 
distribution of a nearby molecular cloud and natu- 
rally lead to the formation of pillar-like structures 
within the cloud. The heads of these pillars may 
experience gravitational collapse and undergo star 
formation. Thompson et al. (2002) detected two 
bright stars in the head of the pillar column III in 
M16 using near-infrared observations. 

We attempt to estimate the column density of 
atomic hydrogen (Ah) per pixel across the WFC3 
field of Figure 1 by using the empirical relation 
(Equation 2) of Seward (1999): 

A H = 1.9 • 10 21 A V atoms cm" 2 mag" 1 (2) 

where Ay is computed from the two-dimensional 
map of E(B — V) g in Figure 1 assuming Ay = 3.1 • 
E(B — V) g . The histogram of the pixel log 10 (AH) 
distribution is shown in Figure 6. The majority 
of the pixels have log 10 (AH) in the range of ap- 
proximately 21.9 to 22.1. The rims of the pillars 
in NGC 3603 are characterized by log 10 (AH) — 
22. This value is probably only a lower limit to 
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the true column density of iVn within the pil- 
lars. For example, Thompson et al. (2002) derived 
log 10 (iVH) = 22.8 inside the dark region of the pil- 
lar column III in M16. 

Since the column density of molecular hydrogen 
(Nh 2 ) is one half of the Nh value (Thompson et 
al. 2002), the pillar rims in NGC3603 likely have 
log 10 (./Vh 2 ) — 21.7, which agrees well with the sim- 
ulation results of Mackey & Lim (2010). With this 
column density at the rims, the central H2 column 
density of the pillar, log 10 (A^H 2 ), may be as high 
as about 22.7 to 23 according to Mackey & Lim 
(2010), which agrees with the the column density 
of MM1 and MM2 (log 10 (7V H2 ) = 22.6 - 23) as 
derived by Niirnberger et al. (2002). 

This value is consistent with the H 2 column 
density of bright-rimmed clouds (BRCs) undergo- 
ing star formation. Urquhart et al. (2009) studied 
a number of BRCs, which are isolated molecular 
clouds located on the edges of evolved H 11 regions. 
Some BRCs show evidence of significant interac- 
tion between their molecular gas with the ionizing 
radiation from their nearby Hn regions. These 
BRCs are thus comparable in their properties to 
the ionized heads of molecular clouds. They have 
been found to host star formation activity, most 
likely triggered by the ionizing front coming from 
the nearby Hn regions. The column density of 
molecular hydrogen in BRCs with triggered star 
formation is 20.9 < log 10 (iVH 2 ) < 22.8 (Urquhart 
et al. 2009). 

The similarity in log 10 (iVH 2 ) between star- 
forming BRCs and the pillars in NGC 3603 sug- 
gests that the heads of MM1 and MM2 may be 
sites of star formation. Indeed Caswell et al. 
(1989) and Dc Pree et al. (1999) found methanol 
and water maser sources in the heads of both MM1 
and MM2; these sources are the typical signatures 
of newborn stars obscured by dusty molecular en- 
velopes. The star formation in the pillars is likely 
triggered by photoionization-induced shocks due 
to the expansion of the H 11 region in which the 
HD 97950 cluster resides. The strong Ha and Pa/3 
emission in heads of MM1 and MM2 is another 
indication that the OB stars in the HD 97950 clus- 
ter are ionizing the pillars. The ionization shock 
fronts were detected at mid-infrared wavelengths 
by Niirnberger & Stanke (2003). 



4. Summary 

In the present paper we extend previous stud- 
ies of the differential reddening around the young 
star cluster HD 97950 in the NGC 3603 giant H 11 
region (e.g., Sung & Bessel 2004). We derive a 
two-dimensional map of the color excess of the 
gas around the cluster from HST/WFC3 images, 
which allows us to measure the Ha/Pa/3 flux ra- 
tio and its decrement due to dust extinction (see 
Calzetti et al. 1996). Our reddening map covers an 
area of 4.9 pc x 4.3 pc with a scale of 0.004 pc per 
pixel (assuming that NGC 3603 is at a distance 
of 6.9 kpc). The median value of the color ex- 
cess within the central cavity (r < 1 pc, where 
the HD 97950 cluster resides) is E(B - V) e = 
1.51 ± 0.04 mag. Going from the cavity to either 
the northern or the southern edge of the field of 
view, E(B — V) g is seen to increase from 1.5 mag 
to 2.2 mag at a distance of 2 pc from the clus- 
ter. From the cluster towards either the eastern 
or western edge of the field of view (at a distance 
of 2 pc), E(B — V) g rises from 1.5 mag to 1.6 mag. 
The average uncertainty of the derived color excess 
is about 0.1 mag. 

We find a shell structure 1.2 pc west of the 
cluster with a mean E(B — V) g of 1.59 mag, 
about 0.08 dex higher than in the central cavity. 
We interpret this shell structure as the surface at 
which the expanding gas shell detected by Clayton 
(1986) in the north-south direction interacts with 
a denser molecular cloud. 

The ionizing radiation emitted by the OB stars 
in the HD 97950 cluster is likely to be responsible 
for the formation of the two molecular pillars MM1 
and MM2 seen 1.2 - 2.5 pc south-west and south- 
east of the cluster, respectively. We use our red- 
dening map to estimate the column density of H2 
in MM1 and MM2. We derive log 10 (iVH 2 ) = 21.7 
in the pillars' rims and up to log 10 (-/VH 2 ) = 23 in 
the pillars' center, in agreement with the earlier 
estimates by Niirnberger & Stanke (2003). Based 
on the velocity gradient detected in earlier studies, 
we argue that the pillars are closer to us than the 
ionizing HD 97950 cluster. 

The strong Ha and Pa/3 emission in the heads of 
the pillars MM1 and MM2 traces the ionization of 
the pillar heads by the OB stars in the HD 97950 
cluster. The pillar heads appear to be undergo- 
ing star formation as indicated by the presence 
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of methanol and water maser sources (Caswell et 
al. 1989; De Pree et al. 1999). Such a star forma- 
tion activity is likely triggered by photoionization- 
induced shocks due to the expansion of the H n 
region surrounding the HD 97950 cluster. 

The two-dimensional map of E(B — V) g derived 
in this work paves the way to deredden individ- 
ual stars in the same field of view by providing 
an upper limit to the true stellar reddening. In- 
dividual reddening corrections in regions suffering 
substantial differential reddening are essential for 
analyses of stellar photometry, e.g., for constrain- 
ing the the age spread in the HD 97950 cluster and 
its surroundings. 

The HST/WFC3 multi- wavelength images were 
processed by Max Mutchler (STScI) and were 
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como Beccari. We thank Jay Gallagher, Igna- 
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discussions. We also thank the referee for use- 
ful comments. X.P. acknowledges support in the 
framework of the Excellence Initiative by the Ger- 
man Research Foundation (DFG) through the Hei- 
delberg Graduate School of Fundamental Physics 
(grant number GSC 129/1). This work was par- 
tially supported by the Sonderforschungsbereich 
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B5) of the DFG. 
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Fig. 1. — Upper panel: The two-dimensional map of the gas color excess, E(B — V) g , around the young 
cluster HD 97950 in the giant Hn region NGC 3603. North is up and east is to the left. The image covers 
an area of 4.9 pc x4.3 pc with a scale of 0.004 pc per pixel at a distance of 6.9 kpc. The inner 120 pixels 
(r < 0.52 pc) are masked out due to the luminous OB stars of the HD 97950 cluster, which are saturated 
in the F555W and F814 W images obtained with HST/WFC3. The color excess is significantly lower along 
the east-west direction (1.5 - 1.6 mag) than in the north-south direction (where it is as high as 2.2 mag on 
average). A shell structure is seen westward of the cluster associated with a local enhancement of the color 
excess. The molecular pillars MM1 and MM2 are visible at 1.2 - 2.5 pc southwest and southeast of the 
HD 97950 cluster. Lower panel: The same map in grey scale. The locations of specific sources (Section 3.1) 
associated with a local E{B — V) g enhancement are indicated. 
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Fig. 3. — The radial dependence of E(B — V) g from the cluster core towards east, west, north, and south. 
The solid line is the median E{B — V) g value, while the upper and lower long-dashed lines are the 84 and 
16 percentile values, respectively, of the E(B — V) g distribution as a function of distance from the cluster. 
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Fig. 4. — Color composite image of NGC3603, obtained with the map of the color excess (red), the Pa(3 
emission image (green), and the Ha emission image (blue). The positions of the shell and of the pillars MM1 
and MM2 are indicated in the image. 
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Fig. 5. — A cartoon to illustrate why the observed E(B — V) g in the rims of the pillar is larger than in the 
center. The pillar is represented by a cylinder (neutral in its interior and with a density gradient from head 
to tail) surrounded by a outer layer of gas ionized by the OB stars of the HD 97950 cluster. See the text for 
a more detailed description. 
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Fig. 6. — Histogram of the pixel column density of atomic hydrogen, which is derived from the color excess 
of each pixel through the empirical relation of Seward (1999). 
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Table 1 

wfc3 exposure times per filter 



Filter 


Exposure time 


F555W 


1000s 


F656N 


1550s 


F814W 


990s 


F127M 


2397.697s 


F128N 


1197.694s 


F139M 


2397.697s 



Table 2 

E(B - V) g OF specific objects in NGC 3603 



Object Name 


Right Ascension 


Declination 


E{B - V) s 




^ h m s j 


|^ o m s j 


[mag] 


Proplyd 1 


11 15 13.13 


-61 15 50.0 


1.6 


Proplyd 3 


11 15 01.41 


-61 14 45.7 


1.9 


Sher 25: northeast lobe 


11 15 09.71 


-61 15 13.2 


2.0 


Sher 25: southwest lobe 


11 15 05.89 


-61 15 25.0 


1.6 



Note. — The designations Proplyd 1 and 3 are adopted from Brandncr 
et al. (2000), while the coordinates of Sher 25 are from Brandner et al. 
(1997a). All coordinates are J2000 coordinates. 
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